Low-frequency conductivity of a non-degenerate 2D electron liquid 
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We study the conductivity of a nondegenerate 2D electron liquid in a quantizing magnetic field for 
frequencies well below the cyclotron frequency. The conductivity is formed by electron transitions in 
which the energy of a photon goes to the interaction energy of the many-electron system, whereas the 
involved momentum is transferred to quenched disorder. The conductivity peak is non-Lorentzian. 
Its shape depends on the relation between the correlation length r^ of the disorder potential and 
the typical amplitude (5f of vibrations of the electrons about their quasi-equilibrium positions in the 
liquid. The width of the peak is determined by the reciprocal time it takes an electron to move over 
Tc (or the magnetic length I, for r^ <l). In turn, this time is determined by vibrational or diffusive 
motion, depending on the ratio rc/5f. We analyze the tail of the conductivity peak for short-range 
disorder. It is formed by multiple collisions with the disorder potential. We also analyze scattering 
by rare negatively charged traps and show that the conductivity spectrum in this case depends on 
both short- and long-time electron dynamics. 

PACS numbers: 73.23.-b, 73.50.-h, 73.40.Hm 
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I. INTRODUCTION 

In recent years much progress has been made toward 
understanding of transport phenomena in strongly in- 
teracting electron systems. The well-known examples 
are the fractional quantum Hall effect (QHE) and 
metal-insulator transition phenomena in low-density two- 
dimensional electron systems (2DES) in semiconductors 
1^. 2DESs are particularly convenient for investigating 
the electron-electron interaction (EEI) , since the electron 
density n can be varied in broad limits. One of the most 
important effects of the EEI is onset of electron corre- 
lations. The extent to which the system is correlated 
depends on the ratio V of the characteristic Coulomb en- 
ergy e^(7rrt)^/^ to the electron kinetic energy i?kin 



(7rn)i/VSkin 



(1) 



(-Ekin is equal to the biggest of the Fermi energy ep and 
ksT). For F exceeding a critical value Tw a 2DES be- 
comes a Wigner crystal. The parameter Tw is numeri- 
cally large. For low temperatures (E'kin = e_F S> fcsT, in 
which case F = r^), Tyy « 37 whereas for a nonde- 
generate 2DES (ef < fcsT) Fvi/ « 130 g. 

For Tw > F > 1, a 2DES is still strongly correlated, 
but it forms an electron liquid. Based on the success 
of the Fermi liquid theory in describing '^He it is often 
assumed that, in the quantum region ep ^ ^bT, the 
electron liquid can be described by a Fermi liquid, too. 
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However, for very large rs, the system may be more com- 
plicated n, |j. 

For large ksT/ep and for Tw > F > 1, a 2DES is a 
strongly correlated non-Fermi liquid. It should display 
a nonstandard behavior. Experimentally such 2DES has 
been investigated in semiconductor heterostructures j|, 
^ , and in much detail for electrons on the surface of liquid 
helium |Q] . A nondegenerate electron liquid has common 
features with an ordinary liquid. It does not have long- 
range translational order and displays self-diffusion, as 
seen in various numerical simulations ||^, ^ . 

In contrast to ordinary liquids driven by forces applied 
to their surface, a 2DES is driven by a "volume" field, 
that is a field experienced by each electron. It comes 
from external disorder, such as defects and phonons in 
semiconductors, or helium vapor atoms and surface cap- 
illary waves on helium. Because the total momentum of 
a 2D electron liquid is changed through "volume" rather 
than "surface" scattering, electron transport is different 
from transport in an ordinary liquid. It is different also 
from transport in a weakly nonideal electron gas. Even 
though the conductivity is of a metallic rather than acti- 
vated type for weak disorder potential, its dependence on 
the parameters of the system, including electron density, 
temperature, frequency, and a magnetic field, should be 
totally different from that for an ideal gas. 

The change of transport coefficients stems from elec- 
tron scattering by disorder being strongly affected by the 
electron-electron interaction in a 2D liquid. During scat- 
tering an electron is coupled to other electrons, and this 
coupling determines the scattering probability. This is 
the case even for a short-range disorder potential, where 
different electrons are scattered off uncorrelated fluctua- 
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tions. In this respect a nondegenerate electron liquid is 
similar to a liquid of vortices in a superconductor. 

In the present paper we investigate the frequency de- 
pendence of the magneto-conductivity axxii-^) of a 2D 
electron liquid. Such dependence is particularly inter- 
esting as it should provide a direct insight into the way 
in which a correlated electron system exchanges momen- 
tum with disorder. Indeed, in the Drude model the low- 
frequency magneto-conductivity transverse to a strong 
magnetic field B is 

dDrl^) — ne^T^^ {uj) /muJl, LOc — eB/rac, (2) 

where lu^ is the cyclotron frequency and is the mo- 
mentum relaxation rate. Here and below we keep in 
axxi'^) the lowest-order term in The frequency dis- 
persion of (TDr comes from the dispersion of Ty^I{uj). In 
the case of single-electron elastic scattering it becomes 
strong for frequencies of order of the duration of a colli- 
sion with a scatterer. 

The single-electron Drude picture does not apply to 
a 2DES if the magnetic field B transverse to the elec- 
tron layer is quantizing. Here, all single-electron states 
except for one or maybe few are localized |11, and 
the single-electron conductivity is equal to zero at zero 
frequency [p^ p^ . Metallic-like conductivity of a nonde- 
generate electron liquid is a result of the electron-electron 
interaction. Earlier we found a way to take this interac- 
tion into account nonperturbatively and calculated the 
static conductivity and the cyclotron resonance for a 
weak short-range disorder potential These results 

were qualitatively and quantitatively confirmed by the 
experiment jl6[ p7| . 

The question of observing the actual dynamics or elec- 
tron scattering in the electron liquid has not been ad- 
dressed previously. From analogy with single-electron 
scattering one may expect that an insight into this dy- 
namics can be gained from the frequency dependence of 
(^xxd^)- Here we develop an appropriate theory and sug- 
gest relevant experiments. We also study the magneto- 
conductivity for two important types of a disorder po- 
tential that have not been discussed for a nondegenerate 
electron liquid: a smooth random potential and a poten- 
tial of rare charged defects. The physics of many-electron 
transport in these cases is significantly different from that 
for short-range disorder. 

A simple argument shows that the frequency dispersion 
of the long- wavelength magneto-conductivity for <C i^c 
is indeed directly related to many-electron effects. Be- 
cause the electron kinetic energy is quantized, the energy 
of an absorbed photon fiuj may go to either the electron 
potential energy in the disorder potential or the Coulomb 
energy of the electron system, or both. For weak disorder 
electrons are not localized, and the disorder potential is 
largely averaged out by electron motion. Then the pho- 
ton energy may only be transferred to the electron sys- 
tem. However, to provide momentum conservation, this 
transfer must be mediated by disorder. 



For short-range disorder, one can think of photon ab- 
sorption as resulting from an electron bouncing off a 
point defect. In a quantizing magnetic field B, a momen- 
tum transfer to a defect 5p leads to an electron displace- 
ment (5r — [c/eB"^) (5p x B. In the presence of radiation, 
this displacement can change the energy of the electron- 
electron interaction by Tito. Therefore by investigating 
the frequency dependence of the absorption cross-section, 
one can find how an electron moves during a collision. 

As we show, for weak short-range disorder the con- 
ductivity Uxxiy) has a peak at w = with a specific 
non-Lorentzian shape, which is fully determined by the 
electron-electron interaction. Of special interest is the 
tail of this peak. For large photon frequencies (but still 
Lo <^ LOc) a large electron displacement 6r is required 
in order to accommodate the energy. Respectively, a 
large momentum has to be transferred to the disorder. It 
may come from multiple electron scattering. The mecha- 
nism has some similarity with that of anomalous diffusion 
transverse to a magnetic field It significantly slows 
down the decay of axx{'-^) with lu compared to the decay 
calculated for single-collision absorption. 

The frequency dependence of axxii-^) is totally different 
for a long-range disorder potential. Of particular inter- 
est is the potential with a correlation length Vc smaller 
than the inter-electron distance. In this case the electron 
scattering cannot be described in the hydrodynamic ap- 
proximation. As we show, here the conductivity displays 
a characteristic cusp at a; = 0. 

In Sec. II we relate the magnetoconductivity of a 
strongly correlated electron liquid to the electron struc- 
ture factor. We introduce the single-site approximation, 
which determines the short-wavelength behavior of the 
structure factor. In Sec. Ill we analyze the frequency de- 
pendence of the magnetoconductivity for a short-range 
disorder potential. We show that, for nonzero frequen- 
cies but <C Wc, the conductivity in quantizing fields 
becomes a nonmonotonic function of B. In Sec. IV and 
the Appendix we develop an appropriate diagrammatic 
technique and study the far frequency tail of the mag- 
netoconductivity due to multiple electron scattering. In 
Sec. V we discuss two types of an intermediate-range dis- 
order: the smooth disorder potential and the potential of 
short-range electron traps; in both cases the conductivity 
is shown to be related to diffusion in the electron liquid. 
Sec. VI contains concluding remarks. 



II. MANY-ELECTRON 
MAGNETOCONDUCTIVITY: GENERAL 
ARGUMENTS 



In the range F ^ 1 the energy of the electron-electron 
interaction (EEI) 



E 



(3) 
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is the largest in the electron system. Therefore even 
where electrons do not form a crystal, F < T\y, electron 
positions r„ are still correlated. The EEI does not change 
the total momentum of the 2DES, and thus does not 
directly affect the long- wavelength conductivity axxi^) 
(the Kohn theorem). However, momentum transfer from 
electrons to defects depends on electron motion, and so 
o'xxi^^) is ultimately determined by the EEI. 

A standard approach to calculating a^x is based on 
finding elementary excitations in the many-electron sys- 
tem and then studying their scattering by a disorder po- 
tential. This approach is not of much help in the case of a 
nondegenerate electron liquid, because elementary exci- 
tations are not known However, for the types of dis- 
order that we are interested in, the frequency-dependent 
conductivity is determined by electron motion on either 
short or long times. This motion can be described even 
when elementary excitations are not known, as explained 
in Appendix [A[ 

We will consider magnetoconductivity in a quantizing 
magnetic field B applied normal to the electron layer, 
exp{hLLjc/ ksT) 3> 1. Then the electron wave function 
is a wave packet. Its typical size is the magnetic length 
I — {h/mtOc)^/'^ . The further analysis is based on the 
observation [ p^ that, in addition to a magnetic field, an 
electron is driven by an electric field Ef from other elec- 
trons. This field is due to electron density fluctuations, 
see Appendix It leads to a semiclassical drift of the 
electron wave packet with a group velocity cEf/B. 

The semiclassical approximation applies for sufficiently 
high temperatures. 



(4) 



Here, is the typical frequency of vibrations of the 
electrons about their quasi-equilibrium positions in the 
electron liquid {ujp is the plasma frequency for B = Q, 
Q. <C i^c)- The picture of moving wave packets, with 
continuous energy spectrum, is qualitatively different 
from the single-electron picture where the electron energy 
spectrum is a set of discrete degenerate Landau levels. 

Electron motion leads to averaging of the disorder po- 
tential. Together with inter-electron energy exchange 
it eliminates single-electron localization by an arbitrar- 
ily weak random potential studied in the QHE theory 
Iprl , [l^ . A typical electron energy in the liquid is fc^T. 
Therefore a sufficiently strong disorder is needed in order 
to localize an appreciable portion of electrons, potentially 
leading to a glass transition. In this paper we assume 
that the disorder potential is weak and the electron liq- 
uid displays self-diffusion and associated self-averaging. 
Specific conditions depend on the correlation length of 
the disorder potential and will be discussed later. 



A. Magnetoconductivity for weak disorder 
potential 

The Hamiltonian of the electron liquid in the presence 
of disorder has the form 



H — Hq + Hee + Hi , 



E 



exp [iqr„ 



(5) 



and Vq are the 



Here, Hq is the sum of the single-particle Hamiltonians 
pl/2m [with p„ = "ihVn + (e/c)A{rn)]; Hee is the 
EEI Hamiltonian and is given by Eq. (| 
Fourier components of the disorder potential 

The long- wavelength magnetoconductivity is given by 
the correlator of the total electron momentum P — J^Pn- 
The latter satisfies the equation of motion 

(e^jy is the antisymmetric permutation tensor). The low- 
frequency conductivity, lu ^ tOc, is determined by slow 
time variation of P. Therefore the time derivative in 
this equation can be ignored. The expression for P can 
be then substituted into the Kubo formula for the con- 
ductivity, giving axxii-^) in terms of the correlator of the 
density operators pq weighted with the disorder poten- 
tial, 



.2,4 



Z4[l-exp(-/3w)] 



OO 



dt 



Aujt 



X^(qq')(VqKj'PqWPq'(0)). (6) 

q.q' 

Here, ( • ) implies thermal averaging followed by averaging 
over realizations of the random potential, S is the area 
of the system, and (3 — fi/ksT. 

In the case of a weak disorder potential, the density- 
density correlator in Eq. (^) can be evaluated to zeroth 
order in Vq (the criteria are discussed below). Then the 
conductivity can be expressed in terms of the dynamical 
structure factor of the electron liquid 



5(q,w) 



dte^'^*5(q,t), 



5(q,i)-A^-i(pq(t)p_q(0))^, 



(7) 



where N — nS is the total number of electrons, and the 
subscript means that the correlator is evaluated in the 
absence of disorder. 

For a nondegenerate liquid, it is convenient to write the 
conductivity (g) in the form of an Einstein-type relation 



ne^Dg 



Ds = fT-\Lj)/4. 



(8) 



Here, Dg can be thought of as a coefficient of electron 
diffusion in the disorder potential; it should not be con- 
fused with the coefficient of self-diffusion in the electron 
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liquid discussed in the Appendix ^). The characteristic 
diffusion length in the disorder potential is given by the 
size of the electron wave packet I, and the collision rate 
is 



q 



where the overlinc denotes averaging over realizations of 
disorder. 

The rate (|^) is quadratic in the disorder poten- 
tial, as in the standard Drude approximation. The ex- 
pression for the conductivity (0), @ goes over into the 
Drude formula (||) if one sets t — Atj^^ / (3uJc- However, 
in contrast to the single-electron Drude approximation, 
the dynamic structure factor in is determined by the 
electron-electron interaction. In particular iS(q, to) de- 
pends on the electron density n. 

The factor [1 — exp{—Puj)]/(3uj in Eq. (||) is equal to 1 
in the most interesting frequency range Plu 1, which 
includes the central part of the peak of the low-frequency 
conductivity. The frequency dependence of the effective 
scattering rate and the conductivity in this range is de- 
termined by S{q,uj). On the other hand, in the analysis 
of the conductivity tail we will be interested primarily 
in the exponent, whereas [1 — exp{—f3uj)]/Puj leads to a 
smooth frequency dependence of the prefactor (oc for 
Poj ^ 1). Therefore we omit this factor in what follows. 



B. The single-site approximation 

The expression is significantly simplified in the im- 
portant and most common situation where the correla- 
tion length of the random potential Tc is small compared 
to the inter-electron distance n~^^^. Here, at most one 
electron is scattered by a given fluctuation of the poten- 
tial, for example, by an impurity in the case of electrons 
in semiconductors or a short-wavelength ripplon in the 
case of electrons on helium. Since the 2DES is strongly 
correlated, all other electrons are far away. 

The condition Tc <C n~^/^ allows us to single out the 
most important terms in the structure factor 5(q, oj). 
The major contribution to the sum over q in Eq. (S) 



comes from q ~ mm{l 



> n^/^. On the other 



hand, S{q,t) as given by Eq. is a double sum of 
exp[iqr^(i)] exp[— ?qr„,(0)] over the electron numbers 
TO, to'. The terms with to 7^ to' are rapidly oscillating 
for q 3> n^^/^. Therefore, when calculating r^^, one 
should keep only diagonal terms with to = to'. We call 
this the single-site approximation. 



5(q,t) 
5ss(q,t) 



( 



.''qr(t)p-jqr(o) 



(10) 



where r = stands for the coordinate of an TOth elec- 
tron. The result is independent of to, and therefore we 



disregarded the electron number in Eq. ([lO|). Respec- 
tively, iS(q, w) ~ 5ss(q, 1^), where 5ss(q, w) is the Fourier 
transform of iSss(q, 

The transition from the Kubo formula (H) to Eqs. (H) 
- ( p^ corresponds to the approximation of independent 
scattering events for each individual electron. It is sim- 
ilar to the standard ladder approximation of the single- 
electron theory and applies provided the duration of a 
collision tcoi is small compared to the reciprocal rate of 
electron scattering by the disorder potential. In turn, icoi 
is the typical time range that contributes to iSss (q, w) . 



III. FREQUENCY DISPERSION OF THE 
CONDUCTIVITY: A SHORT-RANGE 
POTENTIAL 

We will first consider the case of a short-range potential 
with correlation length 



rc « Sf = {kBT/mujlY'^. 



(11) 



Here, 5f is the typical thermal displacement of an elec- 
tron from its quasi-equilibrium position in the electron 
liquid. As explained in Appendix ^ electron motion on 
distances smaller than (5f and for times much smaller than 
is a transverse drift in a nearly uniform fluctuational 
electric field Ef JM [the electron vibration frequency 
is given by Eq. (^]. 

We will calculate the structure factor using a for- 
mulation that differs from the one used in Refs. |l^. It 
is advantageous in that it can be generalized to the case 
of multiple scattering by disorder potential, as shown in 
Sec. IV. 

Thermal averaging of a single-electron operator over 
the states of the many-electron system in Eq. ( p^ ) can 
be done in two steps. First we average over the states of 
a given (TOth) electron for a given many-electron configu- 
ration. Configuration averaging is done next. It comes to 
integration over the relative positions of the guiding cen- 
ters Rm' of all other electrons (to' 7^ to) with respect to 
Rm- The integration has to be done with th e Bol tzmann 
weighting factor exp(— iJee/fcsT) [see Eq. ( All )]. This 
is because the electron kinetic energies are eliminated by 
the Landau quantization, and the only relevant energy of 
the system is the potential energy of the electron-electron 
interaction. 

The first averaging means taking a trace of a corre- 
sponding single-electron operator on the single-electron 
wave functions of an (TOth) electron '0fe(r) (r = r™). The 
functions Vfc belong to the lowest Landau level (LLL) 
and should be found assuming that the electron is in a 
uniform electric field Ef created by other electrons. No 
extra weighting factor (in particular, no Boltzmann fac- 
tor) has to be incorporated when calculating the trace. 
The energy is determined only by the many-electron con- 
figuration, and thermal averaging is done over such con- 
figurations. 
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The wave functions ipk (r) depend on the many-electron 
configurations only in terms of the fluctuational field Ef . 
Therefore the configuration averaging is reduced to aver- 
aging over the distribution of the field Ef (A3), which we 
denote by {•)Et (this notation includes subsequent aver- 
aging over realizations of the random potential, if nec- 
essary). Overall, the average value of a single-electron 
operator 0{r) can be written as 



{0{r)) = 



Y,{Mr)\0{r)\Mr)) 



(12) 



where the prefactor is just the reciprocal number of states 
of the lowest Landau level. 



In order to find the structure factor 5ss(q, t) [Eq. (|10[)] 
to zeroth order in the random potential we will use the 
explicit form of the LLL wave functions ■i/;^,'^-' (r) of an 
electron in the field Ef 



1 



■ exp 



(13) 



Here, we chose the x-axis in the direction of Ef, i.e., 
Ef = E{ x; Ly is the size of the system in the j/-direction, 
and the coordinate x is counted off from —eEi/muj'^. The 
magnetic field is chosen along the negative z-direction, 
B = -Bz. 

The corresponding electron energy e^'^-' (counted off 
from nujc/2 - e2E^/2ma;^)) is 



-(0) 



eEfW. 



(14) 



The structure factor (|T^) is determined by the trace of 
a product of the single-electron operators exp(±jqr(i)) 
taken at different times in the Heisenberg representation. 
From Eqs. (|l|), (|l|) we have 



(0)|p«qr(t)p-»qr(0)|^(0)^ 



k I 



(15) 



where = cEf x B/i?^ is the semiclassical drift ve- 
locity; for chosen axes the vector wd points in the y- 
direction. In Eq. (^5|) we disregarded fast-oscillating 
terms cx exp(±ma;c^) with n > 1. Such terms make ex- 
tremely small contribution to the conductivity for lo <C 



Following the procedure (12), in order to find 5ss(q, t) 
we have to average the right-hand side of Eq. (|l^) over 
the fiuctuational field Ef . For the Gaussian ||l^ field dis- 
tribution ( [A^ ) this gives 



5ss(q,i) = cxp(-ig2/2 



1 



(16) 



Eq. (^^ shows that the structure factor decays very fast 
for wave numbers q ^ 1/L For ql ^ \ \\, also rapidly 
decays with time for t ^ t^. 

The characteristic time is simply related to the r.m.s. 
drift velocity and the r.m.s. fluctuational field (E^^ 



f /' 



V2I 



V2h 



e/(E2) 



1/2 



(17) 



(the choice of the coefficients is convenient for Eq. ( |21 
below). A closely related time 



-1/2 



(18) 



was introduced previously [ [l5[ as the average time of 
flight of an electron wave packet over a distance / in the 
crossed flelds Ef , B. In the case of a short-range disor- 
der potential (rc <C I), te is the duration of an electron 
collision with the fluctuation of the potential (a point 
defect). For a Gaussian distribution of the fluctuational 
field assumed here we have ie = {2/ny/'^te. 

Eqs. (||), (|l6|) give the conductivity in a simple form. 
The typical values of q transferred in an electron collision 
are qc = (r^ -I- Z^)~^/^. The duration of a collision is 



the time over which the correlator (16) decays, te/qj- 
For a short-range potential (^ij), we have te/qj ^ fi^^. 
This justifies the assumption that the fiuctuational field 
remains constant during a collision, which is equivalent 
to the assumption that t ^ il^^. 

The condition for the disorder potential to be weak so 
that collisions occur successively in time is 



[(re/0' + l]'/'«r(0) 



(19) 



where t~^{uj) is given by Eq. (^ and is quadratic in the 
potential strength. 

The typical frequency fl [Eq. (^] of electron vibra- 
tions about its quasi-equilibrium position is also the rate 
of inter-electron energy exchange. The time interval be- 
tween successive collisions of a vibrating electron with the 
same short-range scatterer ~ {qc5{)^fl~^ is much larger 
than Therefore an electron looses coherence be- 

tween successive collisions. This shows that interference 
effects leading to weak localization in the single-electron 
approximation are not important, for weak scattering. 



A. The conductivity for a 5-correlated random 
potential 

The expression for the conductivity can be obtained in 
an explicit form in the important case of a (5-correlated 
random potential V{r)V{r') = w^(5(r — r'), or 



l^qP ^ TTn^j^f/2S, h-y = (2/7r)i/2w/L 



(20) 



The parameter 7 cx 1// introduced here is a convenient 
characteristic of the random potential in the problem of 
electrons in a quantizing magnetic field. It gives the 
width of the peak of the density of states p{E) in the 
single-electron approximation [ pT| ; in particular, on the 
tail p{E) cx exp(-4£'V?iS^), see Ref^. From Eq. M), 
the scattering rate for the potential (pO]) is t~^(0) ~ 7^e, 
and the condition for the random potential to be weak 
takes a simple form 
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Collecting Eqs. (|), (|), ^ we obtain a simple 

explicit expression for the frequency-dependent magne- 
toconductivity, 



o\ (w) 



16 

fi- 



eri (uj) 



Ujtf. 



(21) 



In the quasi-static limit of small ojt 
with our previous result 



Eq. ( pl| ) coincides 
The conductivity (Jxx{Q) 
has a form of a single-electron conductivity in a magnetic 
field, with a scattering rate t'"^ j'^t^ quadratic in the 
disorder potential. However, in contrast to r^/ Eq. (||), 
the value of is fully determined by the EEL Through 
the factor it depends on the fluctuational electric field 
that drives an electron during a collision with a scatterer. 
It scales with the density of the electron liquid as n~^/^, 
so that the overall conductivity axx{0) oc n^^^. 

The frequency dependence of the conductivity is de- 
termined by the dimensionless function (Ti(w), which is 
shown in Fig. |l|. It peaks at zero frequency and monoton- 
ically decays with increasing lo. In contrast to the Drude 
conductivity in the absence of a magnetic field, which has 
a Lorentzian peak a^x 1/(1-1- w^r^), the peak of cti is 
strongly non-Lorentzian. The characteristic width t^^ of 
the peak of ai is independent of the disorder potential. 
Its dependence on the electron density, temperature, and 
the magnetic field is of the form f^^ cx n^/'*T^/^i3~^/^. 




FIG. 1: The frequency dependence of the many-electron con- 
ductivity for short-range disorder as given by Eq. (pl|). 

The tail of the peak is exponential in uj. Interest- 
ingly, it has the exact form of the Urbach rule j2^, i.e., 
jlntJil oc LO. [However, for very large wig where mul- 
tiple scattering becomes important, this form is mod- 
ified, see Sec. 



IV.] The shape of the tail can be un- 



derstood by noticing that the conductivity is formed by 
processes in which the energy Uoj of absorbed photons 
goes to the many-electron system. In an individual ab- 
sorption process the involved electron moves by the dis- 
tance SR = h,uj/e\Ef\ along the fluctuational electric field 
Ef II X. The squared matrix element of a dipolar electron 



transition accompanied by a displacement SR = SR x is 
determined by the (squared) overlap integral of the wave 
functions ( p^ ) 



\{^P^°\r)\^p'^' {r + SR))\'' = exp [-{SRf/2l' 



/,(0) 



(22) 



On the other hand, the probability ( |A3| ) to have a fluc- 
tuational field E{ = hu!/{eSR) is 



p oc exp 



\ e^iSRfi^j))- 



By optimizing the product of the two exponentials 
with respect to SR, we obtain that the conductivity is 
oc exp(— wfe), with account taken of the expression ( p^ ) 
for fg. This is in agreement with Eq. (|2l|). 

It is important to check the assumption that the fleld 
Ef is uniform over relevant distances. In a strongly cor- 



related system |VEt| 



e n 



3/2 



Therefore for optimal SR 



and Ef the relative change of the fleld on the distance SR 
is 



\VEi\5R 



v'/^SR 
~E'f 



\SR)^ 



Tuo 



ilte < 1. 



It is interesting that this condition does not impose lim- 
itations on LO. 

Another interesting feature of the many-electron mi- 
crowave conductivity axxi^) is its nonmonotonic depen- 
dence on the magnetic fleld. Since 7 oc 1/Z oc B^^^ [see 
Eq. (|2^)] and (x B^/^ [Eq. (^], the static conductivity 
Ca;a;(0) oc B^/^ is increasing with B for quantizing flelds. 
This happens because, as B increases, the electron wave 
function becomes more localized, thus increasing the ef- 
fective strength of coupling to short-range scatterers. At 
the same time, the electron drift velocity in the fluctu- 
ational fleld decreases, and as a result the characteristic 
collision duration increases with B, leading to the over- 
all scattering rate ^ 7^ie oc B^^^ [^. The increase 
of 

^xx{Q) with increasing B has been conflrmed experi- 
mentally [p^ . 

For nonzero frequencies, axxi^) displays a peak as 
a function of B, see Fig. |[ For comparatively small 
(but still quantizing) flelds we have Lot^ <C 1, and then 
axxi^^) oc te (X B^/^, as in the static limit. On the other 
hand, for B such that ujte ^ 1, the conductivity falls 
down exponentially with increasing Lotf,. The position 
^max of the peak of the conductivity is given by the 
Golden ratio 



i + Vs 



B 



max ^ ^ 



(23) 



We note that the microwave magneto-conductivity dis- 
plays a peak as a function of B in the single-electron 
approximation as well [Q. However, the shape of the 
peak is totally different. In particular, the high-i? decay 
of the single-electron axx{^) is related to localization of 
electron states in the random potential and is described 



7 




FIG. 2: The dependence of the reduced microwave many- 
electron conductivity a-[{uj) — ijtE(J-i{uj) on the magnetic 
field for weak (5-correlated random potential and for ksT S> 
hu!. Both the characteristic scaling magnetic field Buj = 
B/{Luie)^ and the ratio axx{i^)/o''i{(jj) are independent of B, 
see Eqs. (0) and 



by a power law. Many-electron elfects lead to delocaliza- 
tion, and the decay of the conductivity with B becomes 
exponential, cf. Fig. ||. Of course, for very strong B the 
approximation of weak random potential jt^ <C 1 breaks 
down, and then the decay of axxii-^) with B slows down. 

We note also that in the single-electron theory 
c^xxiQ) = due to electron localization 14|. Therefore 
o'xxi'^) as function of frequency has a peak at nonzero fre- 
quency which depends on B. This is in contrast with the 
monotonic decay of axxii-^) in the many-electron theory 
shown in Fig. ^. 



IV. 



CONDUCTIVITY TAIL: THE EFFECT OF 
MULTIPLE SCATTERING 



It follows from the results of Sec. Ill [see Eq. (13)] that 
electron wave functions in crossed electric and magnetic 
fields display Gaussian decay along the electric field. One 
may expect that multiple scattering of an electron in a 
short-range random potential will slow down this decay, 
as it does for localized electrons in the absence of an 
electric field iQ. In turn, a slower spatial decay may 
lead to a slower decay of the many-electron conductivity 
axxi^) with frequency on the tail where ujif, ^ 1. This is 
because it becomes more probable for an electron to shift 
by a larger distance along the fluctuational electric field 
Ef and therefore to absorb a photon with higher energy. 

The approximation of an electron moving in a static 
uniform field Ef is well justified in the frequency range 
ujte ^ 1- Indeed, the many-electron field Ef changes 
over time ~ see Appendix |a|. It remains constant 

over the duration uj~^ of absorption of a photon, because 
<C ^ The argument in favor of the field 

uniformity is based on the fact that the absorption tail 



is formed by large fields Ef. They are experienced by 
electrons that are far away from their quasi-equilibrium 
positions. The larger is the field the larger the distance to 
the quasi-equilibrium position should be. On the other 
hand, this distance is also the scale on which the field 
is spatially nonuniform. We will see that it will largely 
exceed the electron displacement fiuj / eE{ during absorp- 
tion. Both lengths will be assumed much smaller than 
the inter-electron distance n~^^^. 

A large electron displacement requires a large mo- 
mentum transfer to the random potential, q S> n^^^^. 
Therefore it is a good approximation to evaluate the 
correlator in the expression for the conductivity (^ us- 
ing the single-site approximation. This means that the 
product Pq(i)pq'(0) in Eq. (^) should be replaced by 
Emexp[iqr„(t)] exp[iq'r„(0)]. Then 



<7xx{i^) 



2;4 



dt 



AkBTfi^ 

X E(qq') (^1^1' e'^i--(*)e^'i''-(")) , (24) 
q,q' 

where r = [the result is independent of the electron 
number to]. For Tiui ^ ksT we should replace {kBT)~^ 
with [1 — exp(— /3a;)]/?ia;; however, as noted above, it will 
only affect the prefactor in the conductivity. 

The averaging in Eq. (^) can be done following the 
prescription (|l2|), i.e., one first calculates a trace over 
the single-electron wave functions ij^k in a fluctuational 
field and then averages over the field. In contrast to the 
calculation in Sec. Ill, to allow for multiple scattering 
one should use wave functions found with account taken 
of the disorder potential V{y). In our approximation av- 
eraging over realizations of V{y) and Ef is done indepen- 
dently. It turns out to be more convenient to average 
over V{y) first. 

In what follows the random potential is assumed to 
be Gaussian and (5-correlated, with the correlator (|2^). 
We will use the Green function technique. In contrast 
to what was done in the analysis of the tail of the wave 
function [|l^, this technique has to be formulated in the 
frequency domain. We will show that this leads to a 
somewhat unusual set of diagrams. 



A. The projected Green function 

In order to find the low-frequency conductivity it is 
convenient to use the Green function G'e(r,r') "pro- 
jected" on the lowest Landau level. It is constructed 
from the LLL wave functions V'fc (r) of an electron in the 
random potential V{y) and in the electric field Ef, 



Ge(r,r') = MrWkir'Ke - lO - Su) 



(25) 



where Sk is the single-electron energy of the state k. The 
wave functions are linear combinations of the LLL 
wave functions -f/'i"'' [Eq. (p^)] in the absence of disorder. 
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Following the averaging procedure ([l^), we obtain from 
Eq. (H) 



Re 5: 



de I drdr' 



x{V{r)V{r')dr,A. G,(r,r')G,VJr,r'))E^. (26) 

We will consider the random potential ^(r) as a per- 
turbation. The zcroth-order Green function G'f'' can be 
found using the explicit expressions ([l3|), (Q) for the 
wave functions and the energy in the absence of disor- 
der, 



exp(-[2Z2fc - (x + x') - i{y - y')]^^^) 
e-eEfPk-iO 



(27) 



where the a;-axis is chosen along the field Ef as in 
Eq. (P). 

The function 



g(r,r') = 



1 



2nP 



(28) 



in Eq. ( p7[ ) has a simple meaning. It gives the (minus) 
right-hand side of the Schrodinger equation for the pro- 
jected Green function G£(r, r'), that is it replaces the 
(5- function ^(r— r') in the Schrodinger equation for a stan- 
dard Green function. The function g{r,r') is localized in 
a narrow region |r — r'| < 21 and leads to a Gaussian 

fall-off of Gi°^ for large distances. 

The full Green function G^ is determined by the Dyson 
equation. Its solution can be written symbolically as 

G, = G^°^+G^°'>-V-G^°'>+G^°'>-V-G^°'>-V-G^°'>+. . . (29) 

where the central dot implies integration over internal co- 
ordinates, like / driGf\ri_i,ri)V{ri)G^°\ri,ri+i). We 
emphasize that, even though the Green function G^ is 
projected on the LLL, Eq. ( p9| ) contains the full rather 
than the projected disorder potential V^(r). 

A straightforward calculation shows that, to the lowest 
order in V, the conductivity obtained from Eq. (|2^), co- 
incides with the result of Sec. Ill; in this approximation 
the full Green function G^ in Eq. ( p6| ) has to be replaced 
with G, 



(0) 



B. Diagrams for high-frequency conductivity 



According to Eqs. (p7|), (p8[), the Green function 
Gi"-* (r, r') is mostly localized in a narrow region |r — r'| < 
21. As a function of energy, it peaks at the (scaled) mid- 
point in the Ef-direction, where e = eE{{x + x')/2 for 
y — y' . Near the maximum, 

.9(r,r') 



G(°)(r,r') 



(30) 



With this in mind, we now consider the expression (|2£^ 
for the conductivity in terms of the product of the Green 
functions and think of G^ and G* as given by the pertur- 
bation series (|29|) . In the product of the series we need to 



find terms containing Gf \r,r') and [GfK^{r,r')]* with 



e£:f[r('=^(f,f') -r('=)(r,r')] « huj. 

Such terms describe absorption of a photon accompa- 
nied by an electron displacement by SR — tiw/eEf. The 
displacement results from scattering by the random po- 
tential. 

Graphically, the leading-order contribution to conduc- 
tivity can be represented by a sum of the "fat fish" di- 
agrams illustrated in Fig. || in the coordinate represen- 
tation. In this figure, the wavy lines mark the points r 
and r' where the derivatives of the product of the Green 
functions are taken in Eq. (^6|). The lines above and 
below these points represent the series ( p9| ) for G£(r, r') 
and [Ge+?ia;(r, !■')]*! respectively. The crosses ("x") in- 
dicate the factors V in the expansion (|2^) at points 
where the electron is "scattered". Solid lines between 
the crosses denote the Green functions G^"^ (r^, r^+i) and 
[Gg'^^^(r-, r-_|_;^)]* that describe electron propagation be- 
tween collisions. 




FIG. 3: Diagrams for the conductivity (^6|) in the coordinate 
representation. Segments of the solid lines connecting the 
points r = rs and r' = r'^, in the upper and lower halfplanes 
are the Green functions Gi''' (ri, ri+i) and [G^'^^^(ri, r^+i)]*, 
respectively. Crosses ("x") correspond to the Gaussian 5- 
correlated potential V at points r; [subsequent integration 
over Ti is implied]. Dashed lines indicate averaging over V. 
Wav y l ines mark the "external" vertices V^(r) and V^(r') in 
Eq. ( P6|) and imply differentiation of the product of the Green 
functions over r, r' prior to averaging over V. 

An A'' th order diagram has the total number of crosses 
equal to 2(A-|-1) (two crosses come from the external ver- 
tices in Fig. III). We start with the Green functions that 
correspond to the top and bottom segments in Fig. |^. 
Each of them connects points with equal coordinates, 
r^^^ and ri, respectively. Their product is 

G(°)(ri,ri)[G(%Jr^+i,r^+i)]* 
cx (e - zO - eEfXiy^{e + hu) + iO - eEfXj^^-^y\31) 
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Upon integration over e in Eq. (Pq), the expression 
( ^ ) goes over into 27t'^S(Jilu — eEf{xj^_^-^ — Xi)). This is 
the equation of energy conservation in a photon-induced 
transition ri — > r^r+i. Such a transition is exactly the 
process that diagrams in Fig. ^ describe. 

This picture provides a physical insight into the dia- 
grams. It shows that the absorption does not occur in 
a transition r ^ r' between the external points of the 
diagram. Important for the transition is the admixture 
of the wave functions centered away from r, r' (in fact, 
maximally far away, see below). 

The solid lines in Fig. ^ other than the top and bottom 
segments connect spatially-separated points r^, r^+i. The 
leading exponential terms in the dependence of the Green 
functions on the distance between the points is given by 
overlap functions g{ri,ri+i) [see Eq. (pQ)]. The product 
of the (^-functions entering the A^th order diagrams can 
be written as 

N 

9NiUj}) = 5(1-1, i-i) g(rjv+i, r^+i) n |5(i"^,r,+i)|^ 

1=1 

(X exp(-^(r,-r,+i)V2/2). (32) 

i 

with {tj} = ri, . . . , tat+i. We have incorporated into 
Eq. (|3|) the terms from the "external" coordinates r = 
Ts, r' EE r(,,. 

For large hcu/NeEil the integral over r2, . . . , rjv-n '^^^^ 
be evaluated by steepest descent, with the constraint 
hu! = eE{{xpf_^_l — xi) [the integral over ri cancels the 
factor S^^ in Eq. (p6|)]. As in the case of underbarrier 
tunneling in a magnetic field , the extreme points are 
equidistant, 



= xi + {i- l)huj/NeEf, 



(e) (e) 

Vi -Vi 



for z = 2, . . . , A^ + 1. Then, at the extremum, 
\n~g^({vf})^-{r^leEilf/2N. 



(33) 



(34) 



31) 



This expression has a simple physical meaning. The fac- 
tor A^ in the denominator shows the scattering-induced 
increase of the overlap integral between the electron wave 
functions centered at points ri and r^y^j^, which are sep- 
arated by huj/eEf. 

Except for the top and bottom Green functions 
all other Green's functions in the diagrams in Fig. y are 
nonresonant. Their saddle-point values can be evaluated 
using the approximate expressions ( |30| ) . The diagrams in 
Fig. ||(a) have Ts = r that differs from r^' = r'. These di- 
agrams describe interference of different tunneling paths 
and can be negative or positive. The total of all diagrams 
in Fig. 1^ is, of course, positive. The interference of paths 
affects only the prefactor in the conductivity. The lead- 
ing exponential dependence on the distance huj/eEf and 
on the diagram order A^ is not affected. It is the same for 
all diagrams in Fig. ^ Therefore here we will give results 
only for the "diagonal" terms, which are described by the 
diagrams in Fig. ||(b). 



The A''th -order diagrams in Fig. ^(b) contain a mul- 
tiplier [ttH^^^P /2]^^^ from the intensity of the random 
potential (|2^) . Together with the factors from the energy 
denominators [see Eq. (|30|)] and the factors from the in- 
tegration over dri around the extreme points (^), these 
coefficients give an A^-dependent factor 



Cs = [N-f/ujY'^ [{2s - 3)!! {2N - 2s + 1)!!] 



(35) 



which also depends on the position s (r = r' = r^) of the 
wavy lines in Fig. ^(b) (we assume that l<s<A^ + l). 
For large A^ ^ 1, the factor Cg is maximal for s ~ N/2, 



lnC„ 



2N\n{-f/uj) + 2N 



(36) 



The condition s — N/2 indicates that photon- induced 
transitions preferably occur between the states that are 
maximally (and equally) separated from the "external" 
points r,r'. This is optimal in terms of maximizing the 
overlap integral of states with given energy separation. 

A more detailed calculation [Q shows that, in order 
to allow for compensation from diagrams in Fig. ^(a), 
it is necessary to incorporate corrections to the leading- 
order steepest descent integrals. However, as we already 
mentioned, these corrections do not affect the leading 
term in Inaxx- We also note that, for a given Ef, the in- 
plane conductivity becomes anisotropic. This anisotropy 
leads to different prefactors for the conductivity in the 
directions parallel and perpendicular to Ef [e3. 



C. Frequency dependence of the logarithm of the 
conductivity 

The logarithm of the conductivity axxi^) is given by 
the maximal value of the sum of the expressions (^), ( |3^ ) 
with respect to the order of the diagram A^. To average 
over the fiuctuational electric field, one has to add the 
logarithm of the field distribution and find the maximum 
of the resulting expression over Ef . With the Gaussian 
field distribution ([A^), we obtain 

Inaxxiu^) ~ -(3/21/3) (cuQ^^ [\n{Lu/^) - 1]'^' , (37) 



where te is given by Eq. 

Eq. ( |37| ) is the central result of this section. It shows 
that the logarithm of the conductivity depends on fre- 
quency as a;^/3 for large lu. This form of decay is a re- 
sult of multiple scattering in the random potential, which 
helps an electron move along the fiuctuational field as it 
absorbs a photon. The crossover to Eq. (37) from the 
single-scattering approximation ( pl| ) [where | In axx \ oc uj] 
occurs when the optimal number of scattering events 
A^ ~ [wte/ ln(tt'/7)]^/3 becomes large. The correspond- 
ing value of CO depends on both the many-electron fiuc- 
tuational field and the intensity of the random potential. 

Several comments need to be made about this result. 
First, the diagrams in Fig. ^ are unusual for a trans- 
port problem. They are neither the standard "ladder di- 
agrams" , nor the maximally crossing diagrams. Rather 
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the diagrams in Fig. ^(b) are the maximally wrapped 
(embedded) diagrams for the self-energy. Here they ap- 
pear, because absorption of a photon is accompanied by 
only one "real" scattering by the random potential. The 
role of multiple scattering is to alleviate the Gaussian 
decay of the electron wave function along the fluctua- 
tional field. Second, as the frequency increases, the prob- 
ability of a realization of an optimal fluctuational field 
E{ may become non-Gaussian. In particular, in the re- 
gion where hip(E{) is nearly linear in E{ |]To|, we obtain 



|ln( 



V. 



i.e., even a slower decay than Eq. (p7| 



MAGNETOCONDUCTIVITY IN A SMOOTH 
RANDOM POTENTIAL 



In many physically interesting cases the correlation 
length of the random potential is large (or effectively 
large, see below) compared to the typical size of the elec- 
tron wave packet. A well-known example is provided 
by electron systems in semiconductor heterostructures, 
where much of the disorder potential comes from the 
donors that are spatially separated from the 2DES. A 
sufficiently weak smooth random potential does not lead 
to a glass transition in an electron liquid. The liquid 
should then display a nonzero static conductivity cr^j. (0) . 

We are interested in the effect on transport of the dy- 
namics of individual electrons in the electron liquid. Re- 
spectively, we will consider scattering with momentum 
transfer that largely exceeds the reciprocal inter-electron 
distance, q ^ n~^/^. Such scattering is usually more im- 
portant for magnetoconductivity. If, on the other hand, 
the random potential is smooth on the scale of inter- 
electron distance, Tc » n~^/^, the magnetoconductivity 
can be analyzed in the magneto-hydrodynamic approx- 
imation by considering long-wavelength hydrodynamic 
modes of a viscous electron liquid in a magnetic field 
and disorder potential. 



A. Gaussian potential with correlation length 



rc < n' 



-1/2 



We start with the case of a weak Gaussian random 
potential with correlation length small compared to the 
inter-electron distance. Here, to the second order in ^(r) 
the conductivity is given by Eqs. (^, (||), with the struc- 
ture factor 5(q, w) evaluated in the single-site approxi- 
mation (^. 

The results become particularly interesting and in- 
structive if the correlation length of the potential satisfies 
the condition 



(38) 



where 6{ is thermal electron displacement from a quasi- 
equilibrium position in the liquid (11). 



For a random potential that satisfies the inequality 
(p8|), the structure factor 5ss(q, i) needs to be calculated 



for q l/5{. In other words, the electron displacement 
\^m{t) — rm(0)| in Eq. ( p^ should exceed 6i. Such dis- 
placements occur on times t that largely exceed the recip- 
rocal frequency of vibrations about a quasi-equilibrium 
position il^^. They are due to self-diffusion in the elec- 
tron liquid, with the diffusion coefficient Dee ( A14 ) in- 
troduced in Appendix. In the diffusion approximation 
we have 



5s.(q,t) «exp(-Dee9'|i|), |t| > 



(39) 



The physical picture of scattering by a smooth ran- 
dom potential V{r) is as follows. The guiding center of 
the electron cyclotron orbit drifts transverse to the sum of 
the many-electron fluctuational force — eEf and —'VV{r). 
The field Ef leads primarily to vibrations about a quasi- 
equilibrium electron position. In turn, these vibrations 
result in partial averaging of the disorder potential. This 
is somewhat similar to motional narrowing in nuclear 
magnetic resonance. The averaging is incomplete because 
of self-diffusion of quasi-equilibrium electron positions. 
Therefore the momentum transferred by the disorder po- 
tential, and thus the conductivity, are determined by the 
diffusion rate. 

Frequency dispersion of the conductivity depends on 
a specific model of the random potential. A model fre- 
quently used in the analysis of scattering of 2DES, in- 
cluding the quantum Hall effect, is a random potential 
with a Gaussian correlator p^ . 



exp{-q\l/2), 



(40) 



V{r)V{r') = {vl/2TTri) exp 



r')V2r-e^ 



From Eqs. (^), (|9|), (|39|), and (|4C|) we obtain for the 
magnetoconductivity 



ne V} 



■cre(w). 



(41) 



cTci'^) = 1 + cos(a)) ( Si{uj) ^ '^j ^ sin((I;)Ci((I' 



where uj — ijjr'^/2Dee, and Si(z) and Ci(z) are the sine 
and cosine integral functions, respectively Eq. ( |4l| ) 
is written for a; > 0; the function UxxWA is even in w. 

The function crQ{uj) is shown in Fig. ]^ It decays mono- 
tonically with the increasing oj and has a Lorentzian-type 
tail, CTg, 2 1 up' for (2; ^ 1. An unusual feature of the 
conductivity is that Gq is linear in \uj\ for small It 
has the form UQiS^) — const x |a;|. This nonanalytic be- 
havior is a consequence of the slow decay of the structure 
factor 5ss(q, t) (By) with time for small wave numbers q. 
We note that Eq. (|l]) does not apply for very small |a;|. 
Indeed, the single-site approximation ( |39| ) is valid only 
for q ^ v}!'^ . This means that Eq. (pi] ) can be used only 
for ^ DeeTi Or a) ^ nr^. 

The typical duration of a collision with the random po- 
tential is given by the time r^/Z?ee it takes for an electron 
to diffuse over the length Tc . The necessary condition for 
applicability of the above theory is that this time be small 
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FIG. 4: Frequency dependence of the reduced microwave 
many-electron conductivity ctq^ui) [Eq. dZl])]. The scaled fre- 
quency is ill = Lor^/2Dea. The dashed line shows a Lorentzian 
curve with same maximal value and same asymptotic behav- 
ior for large uj. 



compared to the relaxation time t(0). It sets an upper 
bound on the strength of the disorder potential 



\vg\ < fiDee/l. 



(42) 



The inequality (^2|) can be understood by notic- 
ing that, as mentioned above, a smooth disorder po- 
tential adds to the electron drift velocity a term 
?(c/ei3) V^(z; x q)exp(iqr). Eq. ( ^2| ) corresponds to 
the condition that the root mean square electron dis- 
placement due to this velocity over the collision duration 

/ Dee be small compared to the size of the electron wave 
packet I. Interestingly, the correlation length drops out 
of Eq. @. 

The parameter range where a 2D electron system forms 
a liquid is not limited to temperatures where electron 
motion is semiclassical, ksT 3> Ti^ (^). The results of 
this subsection are valid even for ksT < fifl as long as the 
electron system displays self-diffusion. The parameter 
Si in the inequality (pq) is in this case determined by 
quantum zero-point fluctuations. 

The formalism of this subsection can be applied also to 
the case of a smooth potential created by charged donors 
separated from the electron layer by a spacer of width 
d < The major contribution to the conductivity 

comes from scattering with momentum transfer ^ Ti/d. 
However, the long- wavelength tail of the Coulomb poten- 
tial leads to logarithmic divergence of the relaxation rate 
T~^(0) calculated using Eqs. (|), (H). The analysis of 
this divergence requires a hydrodynamic approach and 
will be carried out in a separate paper. 



B. Electron traps 

An important type of disorder, particularly for MOS 
systems, are electron traps. One can think of them as 
deep short-range potential wells 14r (r) located at random 



positions p^. The overall random potential then is 

nr) = ^Vtr(r-pK) (43) 

The potentials Vtr are not weak, but the trap density 
A^tr is assumed small. In particular, we assume that the 

— 1/2 

inter-trap distance N^_^ ' largely exceeds the correlation 
length of the electron liquid. Then, even though the traps 
will capture some electrons, other electrons will be free to 
move and there still will be self-diffusion in the electron 
liquid. 

To lowest order in iVtr, the magnetoconductivity is 
given by Eq. (|^). For Tiui <C fc^T we write it in the 
form 



2/4 



EE 



n.n' K.K' 



dtf 



x(V„14r[r„(t) - p.] • V„.Ftr[r;(0) - p«-]). (44) 

Because electrons are strongly correlated, only one 
electron may experience the potential of a given trap at 
a time. We assume that this electron gets localized on 
the trap. It then creates a repulsive Coulomb potential 
for other electrons, and they stay away from the trap. 

The picture of one localized electron per trap allows us 
to rewrite the sum of the correlation functions in Eq. ( |4^ ) 
as 

E E(- • •) ^ E^'^'^^t^ • ^'^^^^ [^-(0)])- (45) 

Here, is the coordinate of an electron on a Kth trap 
counted off from the trap position p^. We also disre- 
garded terms with k' ^ k, because different traps are far 
from each other and electrons on different traps are not 
correlated. 

We will assume that a localized electron occupies only 
the ground bound state \g) in the potential Vtr(r) and 
that the energy spacing between the ground and near- 
est excited state is hcog^e 3> ksT. Then electrons very 
rarely escape from the traps or are thermally excited to 
higher states. It should be noted that, in fact, the true 
binding potential is stronger than the "bare" potential 
Vtr, because surrounding electrons contribute to electron 
localization by providing a "caging" potential. The lo- 
calization length of the state \g) in a quantizing magnetic 
field is of the order of the magnetic length I <^ 6i. 

For TiLjg^e S> ksT, thermal averaging of a single- 



electron operator O = 0(i"k) in Eq. (|4^) is done in two 
steps. First, one has to find the diagonal matrix ele- 
ment {g\0\g) on the wave functions of the ground state 
of the trapped electron. Then the matrix element has to 
be averaged over the states of the many-electron system. 
In addition, if we are interested in the conductivity at 
frequencies lo <^ ujg~,e, we have 

{g\V^Vtr [r^m ■ V^^tr [r«(0)] \g) 
^ {g\V^Vt, [r,{t)] \g) {g\V^Vtr [r.(0)] 
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In order to calculate the matrix elements 
(5I V„Vtr(rK)|ff) we note that, besides the potential 
of the defect, a trapped electron experiences a po- 
tential from other electrons. The fluctuating part of 
this potential varies on time n~^. In the range 
of interest, (3V, ^ 1 (@) and f3ujg^e ^ 1, we have 
il <C ujg^e- Therefore a localized electron follows the 
many-electron field adiabatically. The overall force on 
a Kth localized electron is — VVtr (!"«) — eE^;, where E„ 
is the many-electron fluctuational field on this electron 
(see Appendix ^). The diagonal matrix element of 
the force should be equal to zero. This gives for the 
correlator ( |45| ) 

(V.^tr [r^m ■ V.Ftr [r«(0)]) - e^E^{t)E^{Q)). (46) 

For a trapped electron, the behavior of (E„(t)E„(0)) 
is determined by motion of neighboring electrons. As in 
the absence of trapping, this motion is vibrations about 
quasi-equilibrium positions superimposed on diffusion of 
these positions. Because the trapped electron itself does 
not move, the vibrations and the diffusion are somewhat 
different from those in the free electron liquid. However, 
we do not expect that this difference to be too large. 
Indeed, a neighbor of a trapped electron has at most 
one of its six nearest neighbors localized. Therefore the 
correlator (E„(i)EK(0)) should still decay on times of the 
order of the reciprocal vibration frequency 

It is seen from Eqs. (|4|), (^) that axxiO) is deter- 
mined by the integral of the field correlator over time, 
with no extra time-dependent weight. In the semiclassi- 
cal approximation, this integral is simply related to th e 
s elf-d iffusion coefficient Dee of the electron liquid (A12), 
( [Al3| ). For a trapped electron, it should be the same, to 
the order of magnitude, and may only differ by a factor 
^ 1 that depends on F. Therefore 



(0) ~ Nt.e^Dee/ksT. 



(47) 



In contrast to the cases discussed before, the conduc- 
tivity CTa;a;(a;) does not peak at lu — for low frequencies. 
The field power spectrum, which is given by the Fourier 
transform of the correlator (p6|), increases with uj. A 
simple calculation shows that, for a Wigner crystal, this 
increase is linear for uj <ti V. Therefore for an electron 
liquid it should be linear in the range 5f/Dee <^ oj ^ fl. 
The conductivity peaks at ~ 17. The overall width of 
the peak of the low-frequency conductivity is ^ SI. 

Eq. (47) shows that, for negatively charged defects, 
magnetoconductivity is proportional to the defect den- 
sity, rather than the electron density n. For w = 0, it de- 
pends on n only in terms of the self-diffusion coefficient 
Dee- The shape of the peak also depends on n. Therefore 
measurements of the conductivity spectrum should pro- 
vide an insight into both long- and short-time electron 
dynamics in the liquid, i.e., self-diffusion and vibrations 
about quasi-equilibrium positions. We note the similar- 
ity between this problem and the problem of dissipative 
conductivity of a 2D superconducting film with vortices 
in the presence of pinning centers. 



VI. CONCLUSIONS 

In this paper we have found the frequency dependence 
of the conductivity of a nondegenerate electron liquid in 
a quantizing magnetic field for u ujc- We have shown 
that this dependence is extremely sensitive to both short- 
and long-time electron dynamics in the liquid and the 
characteristics of the random potential. 

For a short-range potential, the conductivity is deter- 
mined by large-q electron scattering. It occurs as an elec- 
tron drifts transverse to the magnetic field and the field 
Ef , which is created by density fiuctuations in the liquid. 
The results become particularly simple if the correlation 
length of the potential Tc is less than the magnetic length 
I. Here, the shape of the peak of axx^i-^) depends on one 
dimensionless parameter cute and is given by an explicit 
expression (^). The time te is the time of fiight over 
the length I in the field Ef. It is smaller than the rate of 
inter-electron momentum exchange and depends on 
the electron density, temperature, and magnetic field as 
te cx n' 

-3/4^-1/2^1/2^ Therefore by studying the shape 
of the magnetoconductivity peak for short-range disorder 
one can investigate the short-time electron dynamics as 
a function of the parameters of the electron liquid. 

The tail of the conductivity in the range t^^ <C a; ^ o^c 
is exponential and obeys the Urbach rule. The photon 
energy huj is transferred to the many-electron system via 
a radiation-induced electron displacement along the field 
Ef. The momentum needed for the displacement comes 
from scattering by a fluctuation of the disorder potential. 
As Lute increases, it becomes more probable for an elec- 
tron to experience multiple scattering. This leads to the 
change of the asymptotic behavior from | In axx | cx a; to 
|ln(T2;a;| cx w^/^. Such behavior is described by unusual 
"fat fish" diagrams, which correspond to maximally em- 
bedded diagrams in the self-energy (see Fig. 

The conductivity has a different form for a long-range 
disorder potential. Of particular interest is the case 
where the correlation length Tc exceeds the root mean 
square electron displacement from a quasi-equilibrium 
position in the liquid Si. Here, axxi^) is determined by 
electron diffusion in the liquid. It has a simple form given 
by Eq. (|4l|). The shape of the spectral peak of axx{^) 
depends on one dimensionless parameter ur'^/Dee- It 
displays a smeared cusp at w = and decays as w"^ for 
large ujrH Dee- 
Yet another behavior arises in the case where scatterers 
are short-range electron traps. If the density of trapped 
electrons is small, the 2DES remains a liquid. The con- 
ductivity is expressed in terms of the power spectrum 
of the ffuctuational field Ef that the liquid exerts on a 
trapped electron. The spectral peak of the conductivity 
is located at a frequency ^ Vl. The low- frequency part of 
the peak is determined by electron diffusion, whereas the 
shape of the peak near its maximum depends on vibra- 
tions of electrons about their quasi-equilibrium positions 
in the liquid. 
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APPENDIX A: DYNAMICS OF A 
NONDEGENERATE 2D ELECTRON LIQUID 

A snapshot of a correlated 2D electron liquid is shown 
schematically in Fig. la. In such a liquid, for most of the 
time, the electrons perform small-amplitude vibrations 
about their quasi-equilibrium positions in the potential 
formed by other electrons and the neutralizing back- 
ground. In the absence of a magnetic field, the charac- 
teristic frequency of such vibrations is determined by the 
second derivative of the Coulomb potential at the mean 
inter-electron distance and is given by the short- 



wavelength plasma frequency ujp = {2tt ^ /ra) 



which is the characteristic Debye frequency of a 2D 
Wigner crystal. 



(a) 



co(k) 



(b) 



FIG. 5: (a) A snapshot of a correlated electron fluid (schemat- 
ically). The open circle shows an equilibrium position of one 
of the electrons in the field of other electrons. Because the 
electron is displaced, it experiences a restoring force, which is 
determined by the fluctuational electric field Ef. (b) Phonon 
spectrum of a 2D Wigner crystal in a quantizing magnetic 
field (schematically). 

An important characteristic of electron dynamics is 
the typical fluctuational electron displacement from the 
quasi-equilibrium position (5t in Fig. |^a. If the electron 
motion is classical, i.e., for Tiujp <C fc^T, it can be es- 
timated in the harmonic approximation by setting the 
potential energy e^n^^'^Sf equal to fc^T, 



(Al) 



The necessary condition that 5'^ be much less than the 
squared inter-electron distance is equivalent to F 3> 
1. 

The restoring force on the vibrating electron is given 
by the electric field Ef, see Fig. This field is due 
to electron density fluctuations. In the classical regime, 
these are primarily short-wavelength fluctuations pSl. 



For large F the held should be close to its value estimated 
in the harmonic approximation in electron displacements 
from quasi-equilibrium positions, 



(E?> 



F{V)n^/'^kBT. 



(A2) 



The function F{T) was obtained by Monte Carlo simu- 
lations Ijl^. In the whole range 10 < F < Tw = 130 
it remains essentially constant, varying from 10.5 to 9.1. 
The field distribution p(Et) is Gaussian in the central 
part. 



p(Ef) 



k(E?) 



exp[-E2/(E2 



(A3) 



which is an indication that, in this range of F, electron 
motion is mostly weakly anharmonic vibrations. On the 
far tail, the decay of p(Ef) is slowed down compared to 
Eq. (|A3|) |10 . The fluctuational field is the only charac- 



teristic of the electron liquid, which is needed in order to 
describe the magneto-conductivity a^xi^) in the case of 
weak short-range disorder. 



1. A nondegenerate electron liquid in a strong 
magnetic field 

In a strong magnetic field, 3> cjp, the electron mo- 
tion is separated into cyclotron motion with frequencies 
^ ujc and comparatively slow vibrations of the guiding 
centers about their slowly diffusing quasi-equilibrium po- 
sitions. The coordinates of the guiding centers are 



(A4) 



where, z = — B/_B is the unit vector normal to the elec- 
tron layer, and p is the canonical momentum. From ( ) , 
the components X„ , Y„ obey the commutation relation 



[X„,Y„] = P^h/{mu,). 



(A5) 



The magnetic length I gives the typical size of the electron 
wave packet. 

The dynamics of the guiding centers is described by 
the Hamiltonian of the electron-electron interaction i/ge 
(0) projected on the lowest Landau level. 



Hp, 



(A6) 



Heisenberg equations of motion for R„ can be written 
in a closed form in the important case where the electric 
field on electrons is smooth on the scale I, 



cB"2e„ X B. 



Here, 



E,. 



9R„ 



(A7) 



(A8) 
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is the field on the nth electron created by other electrons 
and calculated ignoring noncommutativity (A_5) of the 
guiding centers' components. 

By linearizing E„ in displacements of the electrons 
from their quasi-cquilibrium positions, one can see that 
the motion of the guiding centers is mostly vibrations 
with typical frequency 



fl = LOp/uJc — 2iTe^n^^'^ /mujc 



(A9) 



This frequency also gives the typical rate of inter-electron 
momentum exchange. In the case where electrons form a 
Wigner crystal, Q is the zone-boundary frequency of the 
lower phonon branch, see Fig. ||b. 

Motion of the guiding centers R„ becomes semiclassi- 
cal for 



(AlO) 



It is determined by thermal fluctuations. From ( AlO| ), 
a typical electron displacement from a quasi-equilibrium 
position S{ given by Eq. ( |Al[ ) is S{ '3> I. Therefore the 
fluctuational electric field Ef, which varies on the dis- 
tance Sf^Js uniform on the magnetic length I, as assumed 
in Eq. (0). 

The probability distribution of the guiding centers 
p(Ri, . . . , Rn, . . .) is given by the Boltzmann equation, 



p(Ri,, 



, Rn 



const X cxp(— iJee/fcsT), (All) 



with Hee given by Eq. (A6). Therefore the results for a 
classical electron liquid in the absence of a magnetic field 
can now be carried over to the case of quantizing mag- 
netic field. In particular, the instantaneous distribution 
of the fluctuational field Ef is the sa me as in the classi- 
cal 2DES for B = 0. The condition (|A10D is much less 



restrictive than hujp <C JtsT. 



because fl <^ uj,j 



For F < 130, Eqs. (|A7|), (|A8|), ( |AllD provide a semi- 
classical description of a nondegenerate electron liquid 
in a strong magnetic field. In the long-wavelength limit, 
this liquid can be alternatively described in the hydro- 
dynamic approximation. The transport coefficients (e.g., 
viscosity) can be found from the correlation functions 
of the liquid (e.g., current-current correlator). In two 
dimensions, because of a large contribution from long- 
wavelength modes, the transport coefficients diverge. A 
self-consistent analysis in the case of a classical 2D liquid 
for small frequencies co and wave numbers q was done 
by Andreev [Eg. He showed that transport coefficients 
diverge as In^'^ lo in the limit of small lu. For a 2D elec- 
tron liquid in a short-wavelength random potential this 
divergence is terminated at the cutoff frequency given by 
the rate of electron scattering in this potential r~^. This 
rate determines the long-wavelength static conductivity 
and thus the decay rate of long-wavelength modes of the 
electron liquid. A corresponding self-consistent analysis 



in the presence of a magnetic field will be given elsewhere. 
A logarithmic correction to the viscosity of a 2D Fermi 
liquid due to electron-impurity scattering was found in 
Ref. ||. 

An important feature of the liquid state, which ulti- 
mately gives rise to a nonzero static conductivity, is self- 
diffusion. The coefficient of self-diffusion D^e can be re- 
lated to the r.m.s. displacement of a particle over a long 
time t > ri~\ 



D,, = Ai?2(i)/4t, ARlit) = {[RJt) - R„(0)]2). 

(A12) 

If we assume that the semiclassical approximation (A7) 
applies for long times t 3> f^^^, then the electron dis- 
placement is related to the correlator of the fluctuational 
electric field by 







ARl{t)^{c/Bf dt' dt"(E„(OE„(i")). (A13) 







The natural scale of the electric field is given by the 
r.m.s. value (A2). Field correlations decay over the time 



Q . Then, the diffusion coefficient becomes 



knT ^ 



miOc 



Dee{T,nt). 



(A14) 



into account that 



Here, we took 

We note that, in the absence of a magnetic field, the 
power spectrum of the field Ef goes to zero for uj ^ 0. In 



this case, the double integral over time in Eq. ( A13 ) gives 
just the average increment of the momentum of an elec- 
tron, which saturates for large t. The electron dynamics 
in a magnetic field is different, and we expect that here 
the integral (A13) linearly increases with t. This con- 
jecture is based on the argument that, because of self- 
diffusion in the electron liquid, AR^{t) should linearly 
increase in time. This increase must be caused by the 
fiuctuational field, since quantum corrections to an elec- 
tron displacement are small [p^ . 

An analysis based on magneto-hydrodynamics with a 
frequency-independent viscosity coefficient leads to an 
extra factor Int in the time dependence of AR'^^{t). A 
similar factor should arise in the function D^e for i <C r. 
However, it should saturate and become a constant for 
t > T. For a correlated liquid, we expect that the factor 
Dee is not large. Dee ^ 1- 

In the absence of a magnetic field, straightforward scal- 
ing arguments give for the self-diffusion coefficient an ex- 
pression of the type ( A14), with ojc replaced by the char- 
acteristic vibration frequency LOp. The numerical factor 
which stands for Dee in this expression is known from 
the data of simulations ^, ^, ^0|, it is ~ 0.1 close to 
the melting transition and increases with temperature. 
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